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Abstract 

Monte-Carlo  simulations  featuring  Kolmogorov  phase  screens  have  long  been  used  in  understanding  both  the  effects 
of  atmospheric  turbulence  on  imaging  and  the  performance  of  turbulence  mitigation  methods.  We  have  developed  a 
new,  high  fidelity,  non-Kolmogorov  phase  screen  model  that  allows  the  orientation  of  anisotropy  to  vary  arbitrarily. 
Our  model  is  used  to  generate  simulated  video  sequences  featuring  a  common  satellite-imaging  target  in  low- 
elevation,  daylight  imaging  scenarios,  where  these  effects  are  thought  to  be  common.  These  sequences  demonstrate 
the  effects  of  anisotropic,  non-Kolmogorov  turbulence,  which  can  be  unintuitive.  We  also  briefly  examine  how 
these  non-traditional  features  affect  the  performance  of  certain  post-processing  techniques,  such  as  speckle  imaging. 


1.  INTRODUCTION 

In  Daylight  Imaging  (DI)  of  Space  Situational  Awareness  (SSA)  targets  many  objects  of  interest  are  often  visible 
from  ground-based  imaging  assets  only  at  low  elevation  angles.  Imaging  at  low  elevation  angles  is  challenging  for  a 
number  of  reasons  and  even  more  so  during  daytime.  By  far  the  most  challenging  aspect  is  the  length  of  the  ground- 
to-space  observation  path  through  the  atmosphere.  Considering  the  long  path  length  turbulence  is  certain  to  be 
severe  confounding  image  recovery  algorithms.  Also  challenging  is  the  uncertainty  regarding  the  nature  of 
turbulence  along  these  paths.  Recent  evidence  [1-3]  suggests  that  turbulence  along  these  paths  does  not  obey 
traditional  Kolmogorov,  5/3  power-law  fluctuation  statistics.  This  work  also  suggests  that  turbulence  may  be  highly 
anisotropic  [3].  The  deviation  of  turbulence  from  these  two  norms  has  been  explored  in  depth  analytically  by  a 
number  of  authors  and  accurate,  high  fidelity  modeling  on  these  effects  via  simulation  is  active  area  of  research. 

In  this  work  we  focus,  primarily,  on  visualizing  the  effects  of  non-Kolmogorov  and  anisotropic  turbulence  on 
imaging  during  DI  engagement  scenarios  using  sequences  of  simulated  images.  The  key  development  enabling  this 
comparison  is  a  high-fidelity  non-Kolmogorov  phase  screen  model.  Monte-Carlo  phase  screen  simulations  are 
indispensable  in  both  understanding  the  effect  of  turbulence,  of  any  kind,  on  imaging  system  and  especially  in 
understanding  the  limits  of  turbulence  mitigation  methods  such  as  Adaptive  Optics,  and  image  reconstruction 
methods  such  as  Speckle  Imaging  (SI)  [4]  and  Multi-Frame  Blind  Deconvolution  (MFBD)  [5].  To  that  point,  we 
briefly  examine  the  effect  of  anisotropy  and  changing  power  law  exponent  on  the  ability  of  SI  methods  to  recover 
images  from  multiple  turbulence  corrupted  images.  In  general,  we  find  that  smaller  power-law  exponent  increase 
image  blur  though  images  are  subject  to  relative  less  tip-tilt  distortion.  On  the  other  hand,  power-law  exponents 
larger  than  the  standard  5/3-value  result  in  increased  tip-tilt  with  less  overall  blurring.  As  one  might  expect,  the 
effect  of  anisotropy  is  to  preferentially  blur  images  along  the  direction  of  the  anisotropy.  Images  recovered  using 
Si-techniques  from  simulated,  turbulence-degraded  imagery  with  these  features  are  surprisingly  similar.  This 
similarity  implies  that  turbulence  strength  and  image  frame  SNR  is  more  critical  than  the  specific  nature  of  the 
turbulence  distortion,  i.e.  power-law  exponent  or  anisotropy. 

2.  BACKGROUND 

Monte-Carlo  simulation  involving  Kolmogorov  phase  screens  have  been  use  for  many  years  and  are  a  standard  tool 
in  atmospheric  optics.  Though  other  methods  exist  [6,7]  most  commercially  available  atmospheric  optics  products  , 
such  as  WaveProp,  make  use  of  phase  screens  generated  using  FFT-based  techniques.  FFT-based  methods  begin 
with  a  matrix,  Wn  ,of  complex  circular  Gaussian  random  numbers  with  the  same  number  of  samples  as  the  desired 
phase  screen.  The  inverse  Fourier  transform  of  the  product  between  an  appropriately  sampled  covariance  matrix, 

h=F-1{wn*mK)T1,2r1  (i) 

results  in  two  independent  phase  screens  with  the  correct  spatial  statistics.  In  Eq.  (1),  O(k)  ,  is  the  Power  Spectral 
Density  (PSD)  model  of  the  turbulence  fluctuation  statistics. 
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A  recognized  weakness  of  this  method  is  an  under  sampling  at  low  spatial  frequencies.  This  weakness  is  especially 
worrisome  considering  that  most  of  the  turbulence  energy  in  a  volume  of  atmosphere  is  found  in  the  outer  scales  and 
therefore  also  the  lower  spatial  frequencies.  Lane  [8]  originally  suggested  what  has  become  known  as  the  “sub¬ 
harmonic”  compensation  method.  Lane’s  method  samples  the  lowest  spatial  frequency  region  of  the  PSD  model 
enhancing  the  accuracy  of  the  model.  Frelich  [9],  and  later  Schmidt  [10],  both  describe  an  implementation  of  the 
sub-harmonic  technique  including  3  sub-harmonics,  though  generalizing  the  technique  to  additional  samples  is  a 
trivial  task. 

Apart  from  some  early,  generalized,  work  by  Martin  [11]  and  Coles  [12]  there  has  been  little  attention  paid  to  the 
non-Kolmgorov  power-law  phase  screen  in  the  area  of  optical  wave  propagation  in  the  atmosphere.  Most  work  in 
the  area,  up  to  now,  has  focused  on  the  development  of  analytical  models.  The  first  non-Kolmogorov  phase  screen 
model  was  reported  by  Gudimetla  [13].  Recently,  Roggemann  [14],  and  separately  Bos  [15]  have  presented  work  on 
simulations  involving  non-Kolmogorov  phase  screens.  In  both  cases,  the  user  may  select  the  power  law  exponent  as 
well  as  inner  and  outer  scale.  The  latter  are  modeled  using  a  generalized  form  of  the  analytical  Tatarskii  spectrum  in 
terms  of  radial  wavenumber  K  and  power-law  exponent  CC  is 


O  (K)  =  A(a)p 


exp  (-k-2/^) 

(k2+K T 


(2) 


In  Eq.(2),  A(a)  -  -2— r(«  - 1)«« 
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r  an^ 


is  a  parameter  that  maintains  consistency  between  structure  function 


and  power  spectrum  models  of  index  of  refraction  fluctuations.  The  parameter  /?  is  a  strength  parameter  equivalent 
to  C\  when  #  =  11/3.  The  variables  K0  =  2 K  /  L0  and  Km  —  c(OC)  /  Z0  parameters  for  the  outer-scale,  L0  ,  and 
inner-scale,  l0  ,  respectively  with  c((X)  being  the  inertial-dissipative  consistency  parameter  [18].  Both  Bos  [15]  and 
Roggemann  [14]  use  include  Lane’s  sub-harmonic  compensation  method.  The  model  used  here  adds  the  ability  to 
vary  both  the  degree  and  direction  of  anisotropy.  Anisotropy  is  introduced  following  the  model  described  by 
Wheelon  [16]  and  later  Gudimetla  [17]  via  affine  transform  of  the  spatial  coordinates. 
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Varying  the  direction  of  the  anisotropy  is  enabled  via  a  second  affine  transform  of  the  covariance-sampling  matrix. 

fx>  =  fxco%6  +  fysmQ 

//  =  -/x  sin  #  +  fy  cos  6  ( 

3.  GENERATED  SIMULATED  IMAGERY  USING  PHASE  SCREEN  MODEL 


Fig.  1:  Pristine  SEAS  AT  Image 


Phase  screens  generated  using  the  technique  described  above  were  used  to  generate  simulated  imagery  of  a  common 
“SEASAT”  (Fig.l)  SSA  target  corrupted  by  atmospheric  turbulence.  Simulation  was  carried  out  by  placing  a  single 
phase  screen  across  D=3.5m  imaging  aperture  using  the  technique  described  by  Roggemann  [4].  When  not 
specified,  turbulence  strength  was  fixed  at  D/r0  =10.  Taking  a  new  draw  of  the  circular  complex  Gaussian  random 
matrix,  Wn  ,  creates  new,  independent  phase  screens  and  subsequently,  new,  statistically  independent  images.  By 
placing  these  images  in  sequence  we  can  create  video  sequences  that  either  demonstrate  the  range  of  fluctuations 
under  static  conditions  or  visualize  the  effect  of  a  changing  variable.  Keystones  for  several  of  these  video  sequences 
are  shown  in  the  next  section. 


4.  VIDEO  SEQUENCES:  EFFECT  OF  VARIED  D/  r0  FOR  VARIOUS  POWER  LAW  EXPONENTS 


Fig.  2:  D/rO  =  10  to  40,  a  =  11/3  (Kolmogorov  turbulence) 


Fig.  3:  D/rO  =  10  to  40,  a  =  3.1  (non-Kolmogorov  turbulence,  near  lower  bound) 


Fig.  4:  D/rO  =  10  to  40,  a  =  3.8  (non-Kolmogorov,  near  upper  bound) 


Fig. 2  show  several  frames  from  a  video  sequence  where  the  turbulence  strength,  in  terms  of  D/r0  is  varied  between 
10  and  40  while  the  power-law  exponent  is  fixed  at  a  =  1 1/3,  or  typical  Kolmogorov  turbulence.  The  full  video 
sequence  begins  with  a  relatively  unperturbed  image  that  gradually  becomes  more  distorted  and  also  shifts  randomly 
around  the  image  frame,  both  characteristics  of  turbulence-induced  imaging  distortions.  Fig. 3  is  similar  to  Fig.2  but 
the  phase  screen  power-law  exponent  has  been  changed  to  a  =  3.1.  We  observe  that  the  tip-tilt/random  motion 
effects  are  nearly  absent  in  this  sequence  while  relative  blurring  is  more  severe.  Finally,  in  Fig.4  the  power-law 
exponent  has  been  changed  to  a  =  3.8.  In  Fig.4  both  blurring  and  tip-tilt  distortions  are  more  severe,  the  latter  being 
more  pronounced  that  in  the  other  two  cases.  We  note  briefly  here  that  the  structure  function  definition  of  turbulence 
is  limited  to  a  range  of  3  <  a  <  5.  However,  the  wave  structure  function,  and  all  parameters  derived  therefrom  are 


limited  to  the  range  3  <  a  <  4.  Our  phase  screen  model  uses  the  Fried  parameter  to  set  turbulence  strength.  Of 
course,  the  Fried  parameter  is  defined  as  the  point  where  the  Kolmogorov  wave  structure  function  is  equal  to  6.88. 


Fig.  5:  Degree  of  anisotropy  varied  from  0  to  3  at  D/rO  =  20 


We  now  focus  on  the  topic  of  anisotropic  turbulence,  the  presence  of  which  is  thought  to  be  commonplace  in  the 
upper  atmosphere.  Accordingly,  any  low-elevation  angle,  earth-to-space  imaging  path  will  traverse  long  section 
sections  of  the  upper  atmosphere  and,  therefore,  likely  subject  to  anisotropic  turbulence.  In  Fig. 5  we  demonstrate  the 
effect  degree  of  anisotropy,  8,  from  0  to  3 .  In  our  model  this  variation  results  in  anisotropy  along  the  vertical 
direction.  As  expected,  in  Fig. 5  we  observe  a  very  clear  increase  in  horizontal  blurring  as  the  anisotropy  factor  is 
increased. 

5.  SIMPLE  SI  IMAGE  RESTORATION 

In  this  section,  we  compare  images  reconstructed  using  a  research  SI  processing  technique  from  sequences  of 
turbulence  imagery.  Our  SI  processing  technique  uses  no  SNR  weighting  or  other  optimizations  and  is  identical  to 
the  bispectrum  method  described  in  [19].  As  in  that  work,  5  sub-planes  of  the  bispectrum  are  averaged  to  determine 
the  image  phase  at  each  spatial  frequency.  Our  objective  here  is  to  determine  if  changes  in  power-law  or  anisotropy 
may  adversely  affect  reconstruction  quality.  In  all  cases,  all  parameters  are  fixed  and  100  input  frames  are  used  for 
each  reconstruction.  To  begin,  additive  and  photon-counting  noise  effects  are  not  included. 


Fig.  6:  100-frame  long-exposure  PSF  at  a=3.66  (11/3)  (left).  100  frame  Si-reconstruction  at  D/rO  =  10  (right) 


Fig.  7:  100-frame  long-exposure  PSF  at  a=3.1  (left).  100  frame  Si-reconstruction  at  D/rO  =  10  (right) 


Fig.  8:  100-frame  long-exposure  PSF  at  a=3.8  (left).  100  frame  Si-reconstruction  at  D/rO  =  10  (right) 

We  consider  the  isotropic,  Kolmogorov  results  shown  in  Fig. 6  (right)  as  the  baseline  for  comparison  in  this  section. 
To  that  point,  we  observe  that  the  PSF  in  both  Fig. 7  and  Fig.8  are  wider  than  the  baseline.  Likewise,  in  Fig.8  the 
PSF  is  wider  and,  by  examining  the  video  sequence  a  considerable  larger  bulk  tip-tilt  component  is  observed. 
Comparing  reconstruction  results  (right)  in  each  figure  we  observe,  in  general,  little  difference  in  reconstruction 
quality  between  images.  That  is  to  say,  there  are  no  features  in  the  baseline  that  cannot  be  observed  in  the  other 
images.  On  a  finer  scale,  both  non-Kolmogorov  cases  result  in  a  degradation  of  image  quality.  In  the  “low”  power- 
law  case  (a  =  3.1)  the  image  is  blurred  compared  to  the  baseline.  It  is  possible  that  additional  tuning  of  the  noise 
parameter  in  the  pseudo- Weiner  inverse  filter  may  further  improve  this  image.  In  contrast,  the  “high”  power  law 
case  (a  =  3.8)  contains  imaging  artifacts  common  in  Si-reconstructions  where  the  object  drifts  beyond  the  image 
frame.  It  is  often  possible  to  correct  these  artifacts  so,  while  undesirable,  they  don’t  rule  out  the  use  of  speckle 
imaging  here. 


Fig.  9:  100-frame,  long-exposure  PSF  (left)  and  Si-reconstruction  (right)  for  non-Kolmogorov  turbulence 

with  anisotropy  factor  8=2. 


Fig.  10:  100-frame,  long-exposure  PSF  (left)  and  Si-reconstruction  (right)  for  non-Kolmogorov  turbulence 
with  anisotropy  factor  8=2.  The  direction  of  anisotropy  varies  continuously  from  0  to  180  degrees. 


In  Fig. 9  we  show  an  example  turbulence-degraded  image,  PSF,  and  the  100-frame  reconstruction  for  turbulence 
with  an  anisotropy  factor  of  £  =  2.  In  contrast  to  Fig. 5,  atmospheric  distortions  are  more  severe  in  the  horizontal, 
rather  than  vertical  direction  in  Fig. 9.  Relative  turbulence  strength  is  kept  fixed  at  D/rO  =10.  The  reconstruction  in 
Fig. 9  is  subjectively  equivalent  to,  if  not  slightly  better  than,  the  baseline  reconstruction  in  Fig.5.  This  result  is 
expected  considering  the  use  of  the  affine  transform  in  our  phase  screen  model.  Consequently,  when  anisotropy  is 
increased  along  one  axis  it  is  decreased  in  the  other  effectively  increasing  coherence.  Thus,  we  might  expect 
improved  resolution  in  the  vertical  direction  relative  to  the  isotropic  case.  In  Fig. 9  it  is  difficult  to  conclude  one  way 
or  the  other  if  such  improvement  is  present.  Considering  the  cursory  nature  of  this  exercise  we  will  leave  such 
exploration  for  another  time. 

In  Fig.  10,  like  Fig. 9,  input  images  are  blurred  by  anisotropic  turbulence  with  a  factor  of  £=2.  However,  rather  than 
being  fixed  along  one  axis,  the  orientation  of  the  anisotropy  changes.  Anisoptropy  is  directed  along  the  horizontal 
axis  at  the  start  of  the  sequence,  as  the  sequence  progresses  the  orientation  changes;  falling  along  the  vertical  axis  at 
the  midpoint  and  returning  to  the  horizontal  axis  again  by  the  last  frame.  By  way  of  implementation  the  angle  of 
anisotropy  varies  from  0  to  180  degrees.  Assuming  fixed  anisotropy  one  can  approximate  this  scenario  as  one  where 
the  orientation  of  the  satellite  varies  relative  to  the  anisotropy  over  the  course  an  observation  pass.  Examining  the 
reconstructed  frame  in  Fig.  10  we  observe  both  more  contrast  than  Fig.5  or  Fig. 9  but  also  more  prominent  artifacts. 
We  speculate  that  additional  tuning  of  the  noise  parameter  may  ameliorate  the  differences  producing  an  imaging 
with  both  improved  contrast  and/or  detail  and  less  prominent  processing  artifacts.  Finally,  we  observe  that  the  PSF 
in  Fig.  10  (left)  is  no  smaller  than  the  isotropic,  Kolmogorov  case. 


Fig.  ll:Sample  image  frame  (left)  and  Si-reconstruction  (right)  in  including  Poisson  and  Gaussian  noise  for 

isotropic,  Kolmogorov  turbulence. 


Fig.  12:  Sample  image  frame  (left)  and  Si-reconstruction  (right)  in  including  Poisson  and  Gaussian  noise  for 
anisotropic  turbulence  with  Kolmogorov  power-law  exponent. 


Fig.  11  and  Fig.  12  repeat  the  cases  in  Fig.5  and  Fig.  10  but  include  both  additive  Gaussian  noise  and  Poisson  noise 
associated  with  photo-detection  events.  Additive  noise  in  this  case  models  both  the  background  noise  and  read- 
noise  associated  with  the  detector.  Noise  in  this  case  was  zero  mean  with  a  deviation  of  (Tn=100.  Poisson  noise  was 
modeled  by  using  normalizing  intensity  to  unit  energy  and  scaling  by  K=5,000,000  photo-events,  D/r0=10  in  all 
cases.  Images  were  zero  thresholded  to  approximate  the  decrease  in  SNR  associated  with  increased  background 
illumination  during  DI  scenarios.  In  Fig.l  1  and  Fig.  12  the  image  on  the  left  is  an  example  input  image  and  the  right 


the  SI  output.  Though  a  point  result,  the  anisotropic  results  in  Fig.  12  are  subjectively  better  than  the  isotropic, 
Kolmogorov  case  in  Fig.l  1.  Again,  this  result  is  not  unexpected,  as  anisotropy  results  in  increased  coherence  in  one 
direction.  Fig.  12  also  includes  faint  artifacts  hinting  that  severe  anisotropy  may  pose  problems  for  speckle  imaging. 
In  Fig.  13  and  14,  example  frames  and  point  reconstructions  are  shown  for  power-laws  of  a=3.1  and  3.8  respectively. 
The  reconstruction  in  Fig.  13  is  degraded  compared  to  Fig.l  1.  The  results  here  make  sense  if  we  compare  the  long- 
exposure  PSF  in  Fig. 6  and  7  and  assume  that  the  effect  of  the  smaller  power  law  is  in  effect  the  same  as  an  increase 
in  turbulence  strength.  In  contrast,  we  judge  Fig.l  1  and  Fig.  14  to  be  roughly  equivalent.  Though  the  PSF  in  Fig.8  is 
also  wider  than  Kolmogorov  case,  as  we  discussed  previously,  the  majority  of  the  distortion  here  is  from  tip-tilt.  The 
bispectrum  is  immune  to  in-frame  tip-tilt  distortion,  and  so,  does  not  degrade  image  quality. 

6.  CONCLUSIONS 

We  have  developed  a  new  and  novel  method  of  producing  turbulence  phase  screens  with  several  features  important 
to  modeling  DI  engagements.  We  have  used  this  model  to  produce  video  sequences  featuring  both  non-Kolmogorov 
turbulence  and  turbulence  with  varied  degree  and  direction  of  anisotropy.  The  improved  modeling  capability 
presented  here  will  allow  for  the  simulation  of  complex  engagement  scenarios  where  parameters  such  as  turbulence 
strength,  power-law,  degree  and  direction  of  anisotropy  all  very  simultaneously.  Preliminary  results,  shown  here, 
indicate  that  basic  SI  techniques  remain  effective  when  turbulence  is  non-Kolmogorov  and  for  reasonable  levels  of 
anisotropy  even  in  the  presence  of  Poisson  and  additive  noise  associated  with  increased  background  illumination. 
The  one  possible  exception  being  scenarios  where  high  power-law  exponents  increase  turbulence-induced  in-frame 
jitter.  Optimistically,  this  jitter  could  be  corrected  using  a  Fast- Steering  Mirror  or  by  dropping  frames  where  the 
object  impinges  on  the  edge  of  the  frame. 


Fig.  13:  Sample  image  frame  (left)  and  Si-reconstruction  (right)  in  including  Poisson  noise  for  anisotropic 
turbulence  with  non-Kolmogorov  power-law  exponent  of  a=3.1  . 


Fig.  14:  Sample  image  frame  (left)  and  Si-reconstruction  (right)  in  including  Poisson  noise  for  anisotropic 
turbulence  with  non-Kolmogorov  power-law  exponent  of  a=3.8. 
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